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Abstract 

We show how to provide suitable gauge invariant prescriptions for the classical spatial 
averages (resp. quantum expectation values) that are needed in the evaluation of clas- 
sical (resp. quantum) backreaction effects. We also present examples illustrating how 
the use of gauge invariant prescriptions can avoid interpretation problems and prevent 
misleading conclusions. 



1 Introduction 



The study of the possible dynamical influence of (small) inhomogeneities on the large-scale 
evolution of a cosmological background has recently attracted considerable interest, from 
both a theoretical and a phenomenological point of view. Such a study is of particular 
relevance for the case of inflationary backgrounds, where macroscopic inhomogeneities nec- 
essarily arise as a consequence of the inflationary amplification of the quantum fluctuations 
of the metric and of the matter sources. In that case, the amplified perturbations may 
become large enough to eventually modify the initial cosmological evolution, triggering in 
this way a true backreaction mechanism (as first pointed out and discussed in [1]). But 
inhomogeneities could also affect in a non-trivial way the present cosmological evolution, as 
suggested by recent interpretations of the dark-energy source as the "backreaction effect" 
of appropriately smoothed-out inhomogeneities [2] (see however [3]). 

In order to take into account the influence of such a backreaction on the large-scale 
geometry what is needed, in particular, is the application of a well defined averaging pro- 
cedure for smoothing-out the perturbed (non-homogeneous) geometric parameters. The 
averaging may be referred to a suitably chosen hypersurface, such as the one comoving with 
the gravitational sources (as in [4]), or one of constant energy density, or one of constant 
curvature, etc: the choice will depend, in general, on the specific problem under study. 
What is also needed is a computation of the non-homogeneous parameters up to the second 
perturbative order, at least when the average of the first-order perturbations is vanishing: 
this is indeed what happens in the case of the quantum fluctuations, whose statistical dis- 
tribution is expected to keep its original stochastic properties even after the inflationary 
amplification. 

The computations of the non-homogeneous corrections to a perturbed cosmological back- 
ground, on the other hand, are affected by a well-known ambiguity due to the possible choice 
of different "gauges" , i.e. of different possible parametrizations of the non-homogeneous de- 
grees of freedom, connected among them by coordinate transformations (see e.g. [5]). The 
problem thus arises of using averaging prescriptions which - quite independently of the 
hypersurface (or otherwise-defined portion of space-time) chosen for the smearing of per- 
turbations - turn out to be independent of the particular parametrization, thus providing 
a result which is gauge invariant, at least to second order (see [6] for a recent attempt 
in this direction, focused on the particular case of the de Sitter geometry). The aim of 
this paper is to propose a covariant average prescription satisfying the above property and 
briefly mention some related applications to primordial cosmological backgrounds. Even 
though our approach is in principle non-perturbative, in this paper, for technical reasons, 
we do not develop it explicitly beyond second-order. A detailed study of the implications 
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of our proposal for estimating backreaction effects in a variety of cosmological setups will 
be presented elsewhere. 

The rest of the paper is organized as follows. In Sect. 2, after stressing the difference 
between coordinate reparametrizations and gauge transformations, we discuss on general 
grounds the gauge transformation properties of four-dimensional volume integrals. In Sect. 
3 we concentrate on the quasi-homogeneous cosmological case, introduce a generally co- 
variant and gauge-invariant averaging prescription for non-homogeneous scalar variables, 
and compute its explicit expression - up to second order in the perturbative expansion of 
all background inhomogeneities - for both spatial averages of classical perturbations and 
expectation values of quantized perturbations. Possible examples and applications of these 
results to the computation of backreaction effects are briefly mentioned in Sect. 4. Our 
conclusions are finally summarized in Sect. 5. 



2 Gauge (non)- invar iance of space-time integrals 

In the context of a smoothed space-time geometry, the computation of the "effective" av- 
eraged parameters relies on the integration of their inhomogeneities over an appropriate 
space-time domain (see e.g. [7]). The behaviour of such averaged parameters under gauge- 
transformations is thus determined by the transformation properties of four-dimensional 
integrals defined over the space-time manifold. In order to discuss this point, let us first 
recall the (slightly subtle) difference between general coordinate transformations (GCT) 
and gauge transformations (GT). 

Let us consider a (typically non-homogeneous) scalar field S(x) defined on a generic 
four-dimensional space-time manifold M4 equipped with a Riemannian metric gw, and 
parameterized by a given set of coordinates x M . Under a GCT: 

x - x = f(x), x = f-\x), (2.1) 

S transforms, by definition, as S(x) — > S(x), where: 

S(x) = S(x). (2.2) 

Under the associated GT (or local field reparametrization) - where, by definition, old and 
new fields are evaluated at the same space-time point x - S transforms as S(x) — > S(x), 
where: 

S(x)=S(f-\x)). (2.3) 

Note that this definition applies not only to infinitesimal, but also to finite coordinate 
transformations, and that the standard GCT relation (2.2) is immediately recovered - using 



2 



Eq. (2.1) - when the above equation is evaluated at the point x. In general, S(x) / S(x), 
namely a scalar field is not invariant under a general GT. 



For the metric tensor the GT associated with Eq. (2.1) is similarly defined as g^ u (x 
9^{x), where: 

dx a dx 13 



dfv df v 



9ap(f 



(2.4) 



The subscript appended to the square brackets denotes that, after computing the derivatives 
of the inverse Jacobian matrix (df /dx) -1 , the variable x has to be replaced by f _1 (x). 
Consequently, for g = det g^ v the GT transformation reads y/—g(x) — ► y/—g(x), where: 



\l~9{x) 



dx 



df 



-gif-'ix)), 



(2.5) 



with \dx/df \ the inverse Jacobian determinant associated with the GCT (2.1) 



Consider now the space-time integral of the scalar S over a four-dimensional region £1 
of M.a, spanned by the coordinates x M , and define: 



F(S,n)= [ d 4 xJ-g(x)S(x). 
Jn(x) v 



(2.6) 



Let us now compute its transformation properties under the associated GT that keeps x 
unchanged. Assuming that the region $7 also stays the same: 



F(S,Q.) -► F(S,Q) 



dTxJ-g(x) S(x), 



(2.7) 



where S and g are given by Eqs. (2.3), (2.5), respectively. We now perform a change of 
integration variable from x to x = f~ l {x) (note that x/x!) to obtain: 

df 



F 



(S,fl)= f 



(/(£)) 



d 4 x 



dx 



-g(f(x))S(f(x)). 



Using Eqs. (2.3), (2.5), we are finally led to 

F(S,Q)= i d A xJ-g(x) S(x) = F(S, h) , 
Jnfx) v 



(2.8) 



(2.9) 



where we have introduced the transformed region Q, defined by Q(x) = Q(f(x)). 

Clearly the gauge-transformed integral (2.9) and the original one (2.6) differ since f2 7^ fi. 
Hence, the expression (2.6) is not gauge invariant if we assume - as we did - that the 
integration region is left unaffected by the GT. 

This general result suggests a direct way to obtain gauge-invariant integrals, simply by 
anchoring the definition of the integration region Q, to the gauge choice so that, under a 
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GT, the region O automatically transforms in such a way as to exactly compensate the final 
variation Q, — Q. This can be done as follows: define the integration region 17 of M4 in 
terms of a window function Wq, so that Eq. (2.6) can be rewritten as 

F(S,n)= [ d A xJ-g{x) S{x) = [ d A xJ-g(x) S(x)W n (x), (2.10) 
Jn(x) v JM4 

and assume that Wn(x), rather than being a function given from the outside, is a scalar 
function built out of the dynamical fields at our disposal (explicit examples of step-like and 
deltadike window functions will be discussed below and in the next section). In this case, 
under a GT: 

W n (x) -► W n (x) = W n (r\x)), (2.11) 

and the integration region is consequently transformed from Q(x) to f2(/ _1 (x)). If we now 
repeat the above procedure (performing a GT and a change of integration variable) , starting 
however from the new definition (2.10), and using (2.11), we arrive at the result 

F(S,n)= [ d 4 xJ-g(x)S(x)W n (x) = f d 4 xJ-g(x) S(x) = F(S,il). (2.12) 
JMa Jn(x) v 

The integral (2.10) is thus automatically gauge invariant as a consequence of its invariance 
under GCTs, with the appropriate "gauge-shifting" of the integration region under the GT. 

Let us conclude this section with an explicit example illustrating how the breaking of 
gauge invariance (and general covariance) is related to the transformations properties of the 
boundary of the integration region. We shall consider a stepdike window function, selecting 
a cylinder dike region Q of M4, with temporal boundaries determined by the two spacedike 
hypersurfaces on which a suitable scalar field A(x) (with timedike gradient d^A) assumes 
the constant values A\ and A2; the region is bounded in space by the coordinate condition 
B{x) < ro, where B{x) is a suitable (positive) function of the coordinates with spacedike 
gradient d^B, and tq is a positive constant. More explicitly, we define 

Wn(x) = e(A(x) - A 1 )9(A 2 - A(x))9(r - B(x)) . (2.13) 

where 6 is the Heaviside step function. 

Using this window in the definition (2.10), and applying a GT, we find that the tem- 
poral boundaries of Q are changed because of the gauge shifting A(x) — > A(f~ 1 (x))). The 
transformation of the spatial boundary is similarly controlled by the transformation prop- 
erties of B(x), If B is a scalar we have B{x) — ► B(f^ 1 (x))), the property (2.11) is then 
satisfied, and the integral is gauge invariant. However, for the cosmological backgrounds 
which are of central interest to this paper, all fields are naturally of quasi-homogeneous 
type, and their gradients are typically time-like. In such a context we cannot covariantly 
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define the spatial boundaries as done for the temporal ones, for lack of appropriate fields at 
our disposal. Using for B(x) a function of the coordinates, which does not change under a 
GT, the integral (2.10) fails to be gauge invariant: repeating the same procedure as before 
we find that the final result (2.12), is replaced in this case by: 

F(S,Q)-F(S,Q) = ( d A xJ-g{x) S(x)AW n (x), (2.14) 

where 

AWn(x) = 9(A(x) - A 1 )9{A 2 - A(x)) [6(r - B(f(x))) - 9(r - B{x))\ . (2.15) 

We see that the breaking of gauge invariance is given by an integral over a ring-shaped 
region lying near the space boundary, B = ro, and having a "thickness" Ar controlled by 
the magnitude of the GT. 

We shall now argue that this breaking term tends to be subdominant for large enough 
spatial volumes. The gauge transformations we are interested in are related to the inho- 
mogeneities themselves, which are assumed to be a small perturbation of the cosmological 
background. We may then expand the coordinate transformation with respect to the first- 
order parameter e M as 

of = x fl + e fl (x) + ... (2.16) 

(our considerations can be extended to second order in the perturbations without altering 
the general conclusions). Inserting the Taylor expansion of Eq. (2.15) in Eq. (2.14) we 
obtain: 

F(S,n)-F(S,n) = 
= - J d A x^-gS(x)9(A{x) - A 1 )9(A 2 - A(x))5 (r„ - B(x)) (2.17) 

and we note that the relevant quantity to be considered is the relative magnitude of the 
breaking term, namely: 

AF _ F(S,fl)-F(S,n) 

— = fWT) • (2 ' 18) 

In order to estimate the above expression we can use coordinates in which A(x) and B(x) 
are identified with the time and radial coordinates t and r. The remaining coordinates will 
be like the two angular coordinates (9 and 4>) on the 2-sphere. We then get: 

^ = J d9d<P^-g(x ) S(x y (xo^J ^ dt J d9dcf> jT° dry/-g{x) S(x)^j 

(2.19) 

where (xo) stands for (t, vq, 9, 0), and e r is the radial gauge parameter of the transformation 
(2.16). 
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Barring cancellations, we see that this quantity is of order e r (xq) / d(ro) , where d(ro) is 
the proper size of the spatial region of integration, and e r is determined by some generic 
perturbation 5X = C e X, where C e denotes the Lie derivative, through: 

6 ™ir (2 - 20) 

This means that e r is related to the perturbation 5X/X through a dimensionful factor 
A (possibly controlled either by the characteristic time scale of the background or by the 
characteristic wavelength of the inhomogeneity) . In any case, we expect AF/F to go to 
zero like the ratio of such A to the proper size of the region considered, and thus to vanish 
in the limit of an infinite spatial volume. 



3 A gauge- invariant prescription for spatial averages 
and quantum expectation values 

We have shown, in the previous section, that four-dimensional integrals can be made exactly 
gauge invariant under general (and finite) coordinate transformations by an appropriate 
definition of the integration domain. We have also briefly discussed how the breaking of 
gauge invariance may be related to the choice of the spatial boundaries of integration. Our 
main task, in the present section, will be the application of the above results to the integrals 
which are of central interest to this paper, namely to the average integrals arising in the 
computation of the cosmological backreaction. 

Depending on the context in which the backreaction is approached, there are two types 
of averaging procedure: spatial (or ensemble) average of classical variables, and (vacuum) 
expectation values of quantized fields. In both cases, the associated integrals (and the results 
of the averaging procedure) are not gauge invariant, in general (see e.g. [8, 9, 10], as well as 
the examples of Sect. 4), and one has to face the problem of transforming and physically 
re-interpreting the results from one gauge to another. In order to solve this problem we are 
led to the question: is it possible to define gauge-invariant (spatial averaging) integrals? 

To answer this question we first observe that expectation values of quantum opera- 
tors can be extensively interpreted (and re-written) as spatial integrals weighted by the 
integration volume V, according to the general prescription 

(...) - W d 3 x (...), (3.1) 
Jv 

where, at least for the cosmological backgrounds we are interested in, the integration volume 
extends to all three-dimensional space. 
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This clearly associates expectation values with classical volume averages. Spatial volume 
averages, on the other hand, can be covariantly obtained from the four-dimensional (and 
possibly gauge-invariant) integrals discussed in the previous section simply by using a delta- 
like window function, so as to confine the integration on the spatial hypersurface chosen to 
perform the averaging procedure. 

Such a hypersurface can be selected, as before, as the one on which a suitable scalar 
field A(x) (with time- like gradient d^A) takes the constant value Aq. We can also specify 
the range of integration across this hypersurface by inserting, into the window function, 
a step-like definition of the spatial boundary: for instance, we can impose on a suitable 
(positive) function of the coordinates B(x) (with space-like gradient d^B) to be bounded 
by the constant positive value r$ > 0. This leads to 

W n (x) = S(A(x) - A )9{r - B(x)), (3.2) 

where the delta function localization may also be seen as the limit A2 — ► A\ = Aq in which 
the "thickness" of the step-like time slicing of Eq. (2.13) is sent to zero. 

Using the above window, we are eventually led to define the averaging prescription for 
a scalar object S as 

_ F(S,Q) _ I M4 d*x^S5(A - A ) 6(r - B) 
W{AoM F(1>n) f M4 d* x ^5(A-A )e(r -B) ■ ^ 

We recall (and stress) that the two integrals appearing in this definition are both either 
gauge invariant or not, depending on the transformation properties of the spatial boundary 
determined by B(x). In particular, they are covariant and gauge invariant if B transforms 
under GT as a scalar, while they are not if the definition of the boundary is independent 
of the gauge, as discussed in the previous section. 

In order to explicitly rewrite the above definition in terms of spatial (three-dimensional) 
integrals we can now exploit the properties of the delta function, and perform the time 
integration. Consider, to this purpose, the change of integration variable from t to t, defined 
by t = h(t, x), where the function h is chosen so as to make the scalar field A homogeneous, 
i.e. 

A(h(t,x),x) =A(t,x) = A { -°\t). (3.4) 
After this change of variable the average (3.3) becomes 

j M4 dtd 3 x^f-g(t, x) S(t, x) 6(A<® (t) - A ) 9(r - B(h(t, x),x)) 

W/{Ao,ro} - / 7— Z ' ( 3 - 5 ) 

J M4 dtd^x^-g(t, x) 6(AW (t) - Ao) 9(r - B(h(t, x),x)) 

where we have used the scalar transformation property S(h(t, x),x) = S(t, x), and we have 
absorbed the Jacobian factor \dt/di\ into the definition of g(t,x), according to Eq. (2.5). 
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The delta-function integration can now be easily performed, and gives 

J E d 3 x^-g(t ,x) S(t ,x)9(r - B(h(t ,x),x)) 
{S){A °' r0} = °k A d3x^W^O(r -B(h(t^xJ) ' (3 ' 6) 

where we have called to the time t when (t) takes the constant values Aq. The suffix T,a 
on the integral is there to recall that we are averaging a scalar object S on a section of the 
three-dimensional hypersurface E^, where the given scalar field A(x) takes the constant 
values Aq. 

The gauge invariance of this expression is determined, as in the four-dimensional case, 
by the transformation properties of the spatial boundary, controlled by B(x). We have 
already stressed that in the case of interest to this paper - i.e. for quasi-homogeneous 
cosmological backgrounds - a natural candidate for the scalar field B(x) with space- like 
gradient is missing. Using for B(x) a function of the coordinates which does not change 
under a GT, and following the discussion of Sect. 2, it is clear that the above integrals will 
be gauge invariant only in the limit of an infinite spatial volume. So, to obtain a gauge- 
invariant average prescription, we must consider the limit in which ro — > oo (or B — > 0), 
and the integration extends over all three-dimensional space. In this limit the step-like 
boundary disappears, and we obtain 

_ fj: Ao d?xy/-g(t ,x) S(t ,x) 
h Au d 3 x^/-g{t ,x) 

This represents our gauge invariant prescription for the average of a scalar object S on 
the hypersurface T,a , for any given scalar field A{x) determining a foliation of space-time 
into spatial hypersurfaces E^ with A = const. We should note, in this generalized gauge- 
invariant definition of the average of S, the presence under the integral not of S but of S, 
i.e. of the variable S transformed to the coordinate frame in which A(x) is homogeneous. 

It should be noted, also, that Eq. (3.7) directly generalizes both classical volume aver- 
ages and quantum expectation values. In the quantum case, in particular, Eq. (3.7) can be 
rewritten as 

/C \ (V-g(to,x) s(t ,x)) 

where the brackets without suffix, on the right-hand side of the equality, refer to the standard 
operatorial prescription (3.1) with integration over all three-dimensional space. We stress 
that the two entries of this ratio are not separately gauge invariant, but the ratio itself - 
equivalent to Eq. (3.7) - is indeed invariant, as can be directly verified by applying the 
above definitions (see also the explicit examples of Sect. 4). 
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For a direct check of gauge invariance, as well as for later applications to quasi homoge- 
neous cosmological backgrounds, we will present here an explicit expansion (up to second 
order) of the generalized average (S)a in terms of conventional averages. The subsequent 
discussion can be applied to both the quantum expectation value (3.8) and its classical 
analogue (3.7), for the average of a quasi-homogeneous scalar variable S(x) whose inhomo- 
geneities can be perturbatively expanded as 

S = S {0) + S {1) + S (2) + . . . (3.9) 

Here is the dominant homogeneous part associated to the unperturbed background, 
contains the first-order non-homogeneous corrections, and so on. A similar expansion 
is assumed to apply also to to other relevant variables g^ v and A. We will also assume that 
the conventional average of all first-order terms is vanishing ((S^) = 0, (g^) = 0, etc 
. . . ) , as is indeed the case of inhomogeneities of quantum origin (like those arising from the 
amplification of the vacuum fluctuations in the context of inflationary backgrounds). 

By expanding Eq. (3.8) up to second order, under the above assumptions, we first obtain 
the result 

(S) Ao = + (S {2) ) + J_ (5 (1) (^)W). (3.10) 

We have then to express the transformed (barred) fields, relative to the frame in which 
A is homogeneous, in terms of the original (unbarred) fields, in a general gauge. To this 
purpose let us first consider the "infinitesimal" coordinate transformation parametrized by 
the first-order, e^, and second-order, e^, generators as [11] 



*" + e f 1} + \ (efofl,^ + e f 2) ) + . . . (3.11) 



2 

where 

= ( e (l)>^(D > £ (2) = (eMS) +4>>) < 3 - 12 ) 

(we have explicitly separated the scalar part from the pure transverse vector part e^, e^ 2 ), 
see e.g. [5]). Let us recall, also, that the associated gauge transformation of the scalar field 
S is, to first order, 

5(1) ^ sW=sW-e° {1) S(°\ (3.13) 

and, to second order, 

5 (2) ^ 5(2) = 5(2)_ e i) 5(l)_^ i)+ale(i) ^^(l) 



1 

+ 2 



e° {1) d t (e° {1) SW) + (V (1) + 8%)) d^S^ _ ^(o)] > (3 . 14) 



where dots denote time derivatives, i.e. S = dtS. Similarly, for the first-order transforma- 
tions of y/— using the tensor transformation of the metric (or expanding Eq. (2.5) up to 
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first order) we have 
(v^) (1) 



(v^) (1) = (v^) {1) + (v^) (0) - V 2 £(1) - d t (in (V=y)W) efoj 

(3.15) 

We can now apply the above general transformations of S 1 and y/— g to the particular 
case of the transformation t — > t, used to make A homogeneous (see Eq. (3.4)). To this 
purpose we first expand to second order the transformation t = h(t, x), and invert it, so as 
to recast it in the general form (3.11), namely: 



t = t + ey u + 



2 ( e (i) S " e (i) 



= (l) ^ 2 l^ 1 )"^ 1 ) + e ( 2 )y 
Inserting this transformation into the (expanded) homogeneity condition, 

A(t,x) =A ( -°\t)+A < U(t,x)+A < V(t,x) =A(°\t), 

we are then able to fix the corresponding transformation parameters as follows: 



+ ... 



4(1) 
= « i(0) ' 



"(2) 



A( 2 ) A^At 1 ) 



AW 



2 ' 



= 0. 



(3.16) 



(3.17) 



(3.18) 



Using these particular values into the general transformation laws (3.13)-(3.15) immedi- 
ately gives us the transformed quantities and (\/—g)^ in terms of A and of 
the unbarred fields S and g. Inserting such results into Eq. (3.10), we can finally rewrite 
the generalized average (or expectation value) of S on the spatial hypersurface at constant 
A = Aq, expanded to second order, in explicit form as 



(S) Ao = S<® + {SW) + 



1 



-<s (1) (V=d 



5(0) 



(V=s)<°> N " vv " ' i(0)(^)(0) 

5(0) 



^ (1) (v^) (1) ) 



i(o) 



ft(ln(^)(°)) 



+ 



+ 



5(0) 



2 {a waU)-IJ^-((aW) 2 ) 



2(A(°)y 



+ 



(0) 



_i_ (( ^(i))2 ) _^ (A (2) ) . 

(^(0))2 U I I A {0)\ I 



(3.19) 



Here all terms are evaluated at t = to, and the brackets to the right of this equality denote 
conventional expectation values or classical averages over all three-dimensional space. We 
clearly see that the above prescription strongly depends on the scalar observable A chosen 
to specify the hypersurface to which the averaging is physically referred. However, for any 
given choice of A, we stress that the expression (3.19) is fully independent (up to second 
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order) of the gauge, i.e. (S)a = (S)a , as can be directly checked by applying to S, A and 
^J— g a general coordinate transformation (3.11) with arbitrary values of the parameters 
(see also the explicit examples of Sect. 4). 

This property of gauge invariance allows to evaluate and compare the average of a scalar 
S(x) on different hypersurfaces, defined by different scalars A(x), by solving the dynamics 
in a single gauge. In other words, unlike what is done till now, we do not have to solve a new 
dynamics every time we change the hypersurface of integration. We should note, instead, 
that the result of the conventional average procedure, i.e. (S) = (S^ + + S^) = 
+ (S^) (corresponding to the first two terms to the right of the equality sign), is not 
gauge invariant. 

Equation (3.19) can be rewritten in an equivalent - but somewhat more convenient - 
form as: 



(S)a 



+ 



5 (0) + (A (2) ) + 
"i(0) 



1 



i(0) 



(v^) {0) 
^(ln(v^)(°) 



<A«(v^) (1) >- 



1 



i(°) 



+ ( i(D A (i) 



+ 



1 A(°) 



-((AWf) 



(3.20) 



where: 



A (o) _ §(0) _ ^^4(0) 
AW ' 



1,2, 



A(°) 



(3.21) 



The above formula allows for an easy consistency check. Suppose that S and A are related 
by an arbitrary function S = S(A). In such a case we should find {S)a = S(Aq). Indeed, 
for S = S(A), inserting in Eq. (3.20) the following (easily derivable) relations, 



A« = 0, 

A (o) = 5 "(o) ( i(o) ) 2 ) 



A (2) = -S" {0) (A^) 2 , 

A (i) = s m a^a^ 



(3.22) 



(a prime denotes the derivative of S with respect to its argument), only the first term on 
the r.h.s. survives. 

Let us conclude this section by noting that the integration measure in Eq. (3.7) contains 
the determinant of the full metric of M.^ and not the determinant 7 = det(7^) of the 
intrinsic metric 7^ = g^ of the hypersurface £4. This represents a further difference 
between our prescription and other averages prescriptions commonly used in the literature, 
using the weight VTT[ f° r the integration over the hypersurfaces determined by the chosen 
space-time foliation (see e.g. [4, 7]). 
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It may be noted that a modified averaging procedure with |7| replacing — g in Eq. (3.7) 
is still compatible with gauge invariance: indeed, such a modification can be obtained from 
a covariant four-dimensional integral simply by replacing S(A — A$), in the window (3.2), 
with the following - more complicated but still covariant - window function: 



In that case, by repeating the same procedure as before, and using the relation g 00 g = 7, 
we end up with a result which can be written exactly as in Eq. (3.7), but with —g replaced 



It should be stressed that even in that case a consistent gauge-invariant prescription 
implicitly requires a coordinate transformation (to the barred frame) of the variables to be 
averaged. In that case, after expressing the barred variables in terms of the unbarred ones, 
the average prescription can be written in explicit second-order form as 



(where, to this order, we can use 7 = det(gij)), replacing Eq. (3.20) for the case of the 
window function (3.23). In this convenient form the same consistency checks discussed for 
(3.20) can be easily repeated. 

The two averaging prescriptions (using, respectively, g or 7 in Eq. (3.7)) are clearly 



inequivalent in all cases where y/—g/ VTtI is non-homogeneous. In the rest of this paper 

we will concentrate our discussion on the - apparently more natural, in a covariant four- 
dimensional context - prescription defined in terms of the full metric g^ u , and corresponding 
to Eqs. (3.19), (3.20). 

4 Examples of expectation values for the quantum 
backreaction 

It is well known that the computation of backreaction effects, due to the presence of small 
inhomogeneities perturbing the large-scale cosmological evolution, is plagued by substantial 
ambiguities arising from the gauge dependence of the perturbative approach and of the 
adopted averaging procedure. For the backreaction of the quantum vacuum fluctuations 
in inflationary backgrounds, in particular, there is a controversial literature (see e.g. [12]- 




(3.23) 



by M- 



(S)a 
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[15]), where different results are obtained on the grounds of computations performed with 
different methods in different gauges. 

The aim of this section is present a few simple explicit examples showing that, differently 
from the conventional procedures used in the literature, our generalized average prescription 
(3.19) always gives gauge-invariant results (up to second perturbative order), thus avoiding 
interpretation problems and preventing ambiguous conclusions. 

We shall work in the context of a spatially flat, FRW background geometry, sourced by 
a single scalar field <fi according to the Einstein equations, and we expand our background 
fields {(f), g^u} up to second order in the non-homogeneous perturbations, without fixing 
any gauge, as follows: 

cf>(t,x) = ^°\t) + <pM(t,x) + ^(t,x), (4.1) 
.oo = -1 - 2«« - 2aV 9t0 = -\ + *«) - \ (Pf + BP) , 



9ij = a 2 



+ 2 



Sa (l " ^ - 2^)) + D i3 (EU + + \ ( x (i) + + ,g)) 



(4.2) 



where Dij = didj — 5jj(V 2 /3), and a = a(t) is the scale factor of the homogeneous FRW 
metric. Here f3^\ tp^, are pure scalar first-order perturbations, B^ and xp 

is a trarplpss flnH transvprsp t. 



are transverse vectors (d l B^ = and d l x^p = 0), is a traceless and transverse ten- 



sor {d l h^ = = h^ % ), and the same notation applies to the case of the second-order 
perturbations. 

Let us start with the (almost trivial, but instructive) example concerning the average 
of the scalar field (f), using as our family of reference hypersurfaces the ones where 4> itself 
takes constant values. Let us then apply our average prescription (3.19) with S(x) = <j>{x) 
and A{x) = <f>(x), on the hypersurface where Aq = <^°)(io)- This is a particular example of 
the general case S = S(A) discussed in Eq. (3.22) and it is clear, already from an intuitive 
point of view, that the average should yield the result: 

(<P(t,x))^ = <A (0) (to). (4.3) 



This is indeed the case, as can be quickly checked by computing the average (3.19) in 
the so-called Uniform Field Gauge (UFG) fixed by the condition <f>^ = = (j)^ 2 \ i.e. the 
gauge where the scalar field is homogeneous. In that case, in fact, all terms on the r.h.s. 
of Eq. (3.19) are identically vanishing except the first one, and the result (4.3) is trivially 
obtained. But the same result is obtained in any other gauge, even if <f> is not homogeneus, 
and 4>( l \ (f>^ are both nonvanishing: putting S = (f> = A, in fact, we find that there is an 
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exact cancellation among all terms on the r.h.s. of Eq. (3.19) except the first one, and we 
are led again to the result (4.3). 

The situation is different if we consider instead the conventional average prescription, 
corresponding to the first two terms on the r.h.s. of Eq. (3.19), namely to 

m,x))=^\t) + {^ 2 \t,x)). (4.4) 

In such a case we can obtain different results in different gauges. In the following we 
shall reserve the notation (. . .)^ for our gauge invariant prescription and (. . .) gi for the 
conventional quantum average in the gauge gi. 

In the particular case of the UFG, for instance, we have 4>^ = and we exactly recover 
the previous result (4.3) for any given model of background evolution. But in other gauges 
the correction to the background value <f>^ can be nonzero. We may consider, for instance, 
the so-called Uniform Curvature Gauge (UCG) (also called off-diagonal gauge [16]), where 
we set to zero the variables if) and E of Eq. (4.2) (see also [17] for a general discussion of 
various possible gauges). In that gauge (f)^ / 0, in general, as also confirmed by previous 
studies of the backreaction in models of chaotic inflation [14] with scalar-field potential 
V = m 2 (j) 2 /2. For such models, the computation of the vacuum expectation value (v.e.v.) 
of the quantum fluctuations of the scalar inflaton field in the UCG (neglecting vector and 
tensor perturbations to first order), expanded in the long- wavelength limit, and to leading 
order in the slow-roll parameter e, gives, according to Eq. (4.4), 



(4>(t,x))ucG — 



1_ ^ 



q> 



(4.5) 



(Mp is the Planck mass, and e = —H/H, where H = a/a). This is clearly different from 
the UFG result (<f>(t, x))ufg = <^°H*)> van d for the same model of chaotic inflation to all 
orders in e. There is nothing wrong, of course, with such a discrepancy, since the physical 
meaning of (4>(t, x)) is different in different gauges (0 is not a gauge-invariant perturbation 
even at first order). 

As a second example we will consider the average of a scalar quantity playing a central 
role in the computation of the inflationary backreaction [13]- [15]: the so-called volume 
expansion O, determined by the covariant divergence of a bundle of comoving curves. For 
a scalar-field-dominated geometry, in particular, we have 



6 = VX, u„ = ; , (4.6) 



so that 6 = 3H for the unperturbed homogeneous geometry. 
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If we include perturbations, up to second order as in Eq. (4.1), (4.2), and we compute 
the average of according to the standard prescription, the result is notoriously gauge 
dependent. Consider, for instance, the previous model of chaotic inflation: including the 
contribution of non- homogeneous quantum fluctuations (to second order), and working in 
the UFG, one finds for the v.e.v. of 6 no correction to the background value to first order 
in the slow-roll parameter e, namely: 



(see also [12, 13] for similar UFG results). Repeating the computation in the UCG one finds 
[14], instead, a nonvanishing backreaction even to first order in e: 



variable for scalar perturbations [18]. 

This difficulty (and possible source of misleading conclusions about the backreaction of 
quantum fluctuations in chaotic models) is absent if we apply our prescription (3.19), which 
gives the same result in any gauge. If we refer, in particular, such a covariant average to 
the convenient hypersurface at constant inflaton (i.e. putting A = (f> as before), we find 
that there is no backreaction to first order in e, 



(in agreement with the UFG result of the conventional v.e.v. prescription). We stress that 
the above result is insensitive to the choice of the gauge but not insensitive - in principle - 
to the choice of the reference hypersurface. 

The two examples presented up to now might suggest that the result of a conventional 
average, when performed in the UFG, should coincide with that of the gauge invariant 
procedure performed with A = <f>. More generally, one might even think that the result 
obtained in a gauge where a given scalar variable is homogeneous should coincide with 
the result of the gauge-invariant average performed on the hypersurface where the same 
variable is constant. This, however, is not in general true. Considering the choice A = <p, 
for instance, we find that the covariant average {S) ( f >o reduces in the UFG to the first three 
terms on the r.h.s. of Eq. (3.19): hence, it differs from the conventional average (S) - 
corresponding to the first two terms on the r.h.s. - by the contribution of (S^\y/—g)^). 
Only when this term is zero the two results coincide. 

For a further illustration of this point we may recall that, in the context of the UFG 
choice, we are free to fix the gauge by imposing another condition on the scalar compo- 
nents of the metric perturbations: in particular, we can set to zero one scalar variable 



(@) UFG = 3H l + 0(e 2 ) 



(4.7) 




Note that, in this gauge, also corresponds to the canonical, gauge-invariant Mukhanov 



(e) 0o = 3H 1 + 0(e 2 ) 



(4.9) 



15 



among (3, -0 and E of Eq. (4.2) (but not a, or the background would become trivially 
homogeneous up to first order, see [19]). These three different possibilities identify three 
possible UFGs, connected among them by a coordinate transformation of the type (3.11) 
with eKs = e%s = (namely, by a transformation leaving the time parameter unchanged). 
Under such a transformation is unchanged, but not and (-*/—#) ^ (see Eqs. (3.13)- 
(3.15)). Hence, going from one UFG choice to another, the conventional average (S) (only 
depending on (S^)) may change, with a corresponding change of the term (S^(y/—g)^), 
controlling the difference with the gauge invariant average (obviously, the sum of the two 
terms, determining the gauge-invariant combination (S)^, is left unchanged). 

It seems appropriate to include here an explicit example illustrating the differences 
among the various UFG choices. We shall report the results of a computation [19] of (0) 
for a cosmological background perturbed around a power-law geometry with scale factor 
a(t) ~ l^l 1 / 3 (solution of the Einstein equations with vanishing scalar potential). Working 
in the UFG fixed by = = [3^ one obtains: 

45 (QW^ren 



(@)uFG/3 = 3H 



1 + 



Ml 



while in the UFG with E^ = o 



(6) 



UFGE 



£(2) 

one obtains: 

3(q« 2 ; 



3H 



REN 



(4.10) 



(4.11) 



Here Q denotes the gauge-invariant Mukhanov variable, and the suffix REN denotes that 
the v.e.v. has been regularized through a suitable adiabatic subtraction. These two results 
not only differ between each other, but also differ from the one obtained with the gauge 
invariant prescription (3.19) with A = </>, which gives 



i 



15 (q« 2 ; 



REN 



Ml 



(4.12) 



in any gauge. 



It follows, in particular, that we cannot try to solve the problem of the gauge depen- 
dence of the backreaction considering the UFG as a privileged gauge, as often suggested 
in the literature. The need for a truly gauge invariant approach is evident from the above 
example but, for a better understanding of its importance, let us conclude the section with 
a few comments stressing the different effects of gauge transformations on different types 
of average prescriptions. We shall compare, in particular, our approach with the one more 
frequently used in the literature. 

According to the approach proposed in this paper the averaging procedure is covariantly 
defined, on the hypersurface where a suitably given scalar field A(x) is constant, in terms 
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of fields (S, g^ v , . . . ) evaluated in the coordinate frame (or gauge) where A is homogeneous 
(the barred frame, following the notations of Sect. 3). Only in this frame - i.e. only in 
terms of the barred fields - our average can be written in the simple classical or quantum 
form given, respectively, by Eq. (3.7) or (3.8). If we change gauge, namely if we want to 
rewrite the average in terms of our fields evaluated in a different coordinate frame where A 
is no longer homogeneous, then our prescription assumes a new, more complicated form if 
expressed through the conventional definition of volume averages (or expectation values), 
see Eq.(3.19). However, such a new form is always covariantly related (by definition) to the 
original prescription given in Eq. (3.7) (or (3.8)), and this definitively explains why, at the 
end of our computations, the same result is obtained in any gauge. 

According to the more conventional approaches appearing in the literature, instead, 
when one changes the gauge - i.e. when one moves from the frame where the scalar field 
A{x) is homogeneous to the frame where, for example, another variable (or a set of not 
necessarily scalar variables) is homogeneous - one still maintains the same definition of 
spatial average (or expectation value) used in the old gauge, and simply expresses it in 
terms of the new fields, i.e. of the fields trasformed in the new frame. If the original 
prescription is covariantly defined, and if one of the variables imposed to be homogeneous 
in the new gauge, say Z(x), is a scalar, then the new average prescription can be interpreted 
as the same as the previous one, performed however on a different hypersurface: the one 
where Z is constant (for example, we can interpret in this way the results in [4]). But if 
there is no homogeneous scalar variable available in the new gauge, and/or the average is 
not covariantly defined, then changing the gauge implies a true modification of the averaging 
prescription (not equivalent to a change of the reference hypersurface). In both cases it is 
clear that, within this conventional approach, different gauges lead to different results for 
the averaging procedure. 



5 Summary 

Let us briefly summarize our main results and conclusions. It is generally acknowledged 
that, in order to extract physically meaningful results from cosmological perturbations, it 
is necessary to eliminate all possible gauge artifacts. At linear order in the perturbations, 
the way to do this is well known: following Bardeen's pioneering work [20], it consists of 
defining gauge invariant combinations of the perturbations themselves. Such gauge invari- 
ant combinations have a clear physical meaning, typically representing the perturbation of 
various quantities on hypersurfaces along which some given variable is constant: they are 
like "relative" fluctuations of different variables. 
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A more complicated problem arises when perturbations are considered to second order. 
This can be necessary for a variety of reasons: i) when the perturbations are too large 
to justify the leading order approximation; ii) when the leading-order result vanishes (e.g. 
for non-gaussianity in simple inflationary model); or, finally, iii) when one is interested in 
evaluating backreaction effects either at the classical or at the quantum level. 

The procedure followed at first order can be generalized to define gauge invariant com- 
binations of the second order perturbations (see, for example, [21]). However, whenever one 
has to deal with averages, like in the issue of backreaction effects discussed in this paper, 
such second-order gauge invariant variables are not of much use. 

For the cosmological backreaction, in particular, a simplification occurs due to the fact 
that we are interested in spatial (or quantum) averages of second-order perturbations at very 
large scales. By carefully studying the gauge-transformation properties of four-dimensional 
integrals (while stressing the distinction and relation between gauge and general coordinate 
transformations) and by then going to the limit of an infinitesimal time-like shell, we have 
been able to propose a general formula for the classical or quantum average of any scalar 
quantity, on hypersurfaces on which another given scalar quantity is constant. Our non- 
trivial proposal, Eq.(3.19), can be shown to be gauge-invariant in the quantum case and in 
the classical one for very large spatial-averaging volumes. 

Our proposal is at variance with the recipes that can be found in the literature (often 
with contrasting results and conclusions) since, in any gauge, it contains important extra 
terms. Neglecting those terms induces, in general, gauge-invariance violations that are of 
the same order as the effect one wishes to estimate. This may well explain some of the 
disagreement that can be found in the current literature. 

Unfortunately, the requirement of gauge invariance alone does not fix uniquely our 
prescription. Our simplest gauge-invariant recipe involves the determinant of the four- 
dimensional metric and comes naturally from taking the thin-time-slice limit of a gauge- 
invariant four-dimensional integral. It is also possible to define another gauge-invariant 
average (using Eq. (3.23)), which, in adapted coordinates, reduces to a more conventional 
integration measure involving the induced metric on the hypersurface on which the average 
is computed. While the choice between one or the other definition should be dictated by 
the physics of the problem, in both cases the proposed gauge-invariant averages differ in a 
non-trivial way from the ones that are currently used. 

In this paper we have only given the general ideas behind our proposal, have checked its 
consistency in a number of cases, and discussed some simple applications. A more detailed 
discussion of its implications in specific issues is deferred to further work. 



18 



Acknowledgements 

It is a pleasure to thank for their warm hospitality and their support the following scientific 
institutes, where part of this work was performed: the Galileo Galilei Institute for Theo- 
retical Physics (Florence), the Theory Unit of the Physics Department of CERN (Geneva), 
and the Department of Physics of the University of Bologna. One of us (GV) is also grateful 
to G. Vilkovisky for useful correspondence. 

References 

[1] N. C. Tsamis and R. P. Woodard, Ann. Phys.(N.Y.) 238 (1995) 1; Ann. Phys.(N.Y.) 
253 (1997) 1; V. F. Mukhanov, L. R. Abramo and R. H. Brandenberger, Phys. Rev. 
Lett. 78 (1997) 1624; L. R. Abramo, R. H. Brandenberger and V. F. Mukhanov, Phys. 
Rev. D 56 (1997) 3248. 

[2] S. Rasanen, JCAP 0402 (2004) 003; E. W. Kolb, S. Matarrese, A Notari and A. 
Riotto, Phys. Rev. D 71(2005) 023524; E. W. Kolb, S. Matarrese and A. Riotto, New. 
J. Phys. 8 (2006) 322. 

[3] A. Paranjape and T. P. Singh, Phys. Rev. Lett. 101 (2008) 181101. 

[4] T. Buchert, Gen. Rel. Grav. 32 (2000) 105; Gen. Rel. Grav. 33 (2001) 1381. 

[5] V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rep. 215 (1992) 
203. 

[6] J. Garriga and T. Tanaka, Phys. Rev. D 77 (2008) 024021. 

[7] T. Buchert and M. Carfora, Class. Quantum Grav. 19 (2002) 6109; Phys. Rev. Lett. 
90 (2003) 031101. 

[8] W. G. Unruh, astro-ph/9802323. 

[9] A. Ishibashi, and R. M. Wald, Class. Quant. Grav. 23 (2006) 235. 

[10] S. B. Giddings, D. Marolf and J. B. Hartle, Phys. Rev. D 74 (2006) 064018. 

[11] M. Bruni, S. Matarrese, S. Mollerach and S. Sonego, Class. Quant. Grav. 14 (1997) 
2585. 

[12] L. R. Abramo, and R. P. Woodard, Phys. Rev. D 65 (2002) 063515. 
[13] G. Geshnizjani and R. H. Brandeberger, Phys. Rev. D 66 (2002) 123507. 



19 



[14] F. Finelli, G. Marozzi, G. P. Vacca, and G. Venturi, Phys. Rev. D 69 (2004) 123508. 

[15] L. R. Abramo and R. P. Woodard, Phys. Rev. D 60 (1999) 044010; Phys. Rev. D 
60 (1999) 044011; Phys. Rev. D 65 (2002) 043507; B. Losic and W. G. Unruh, Phys. 
Rev. D 72 (2005) 123510; G. Geshnizjani and N. Afshordi, JCAP 0501 (2005) 011; G. 
Geshnizjani and R. Brandenberger, JCAP 0504 (2005) 006; R. H. Brandenberger and 
C. S. Lam, arXiv:hep-th/0407048; P. Martineau, R. H. Brandenberger, Phys. Rev. D 
72 (2005) 023507; G. Marozzi, Phys. Rev. D 76 (2007) 043504. 

[16] R. Brustein, M. Gasperini, M. Giovannini, V. Mukhanov and G. Veneziano, Phys. 
Rev. D 51 (1995) 6744. 

[17] M. Gasperini, Elements of String Cosmology (Cambridge University Press, Cambridge 
UK, 2007), Chapter 8. 

[18] V. F. Mukhanov, JETP Lett. 41 (1985) 493; Sov. Phys. JETP 67 (1988) 1297. 

[19] M. Gasperini, G. Marozzi and G. Veneziano, in preparation. 

[20] J. M Bardeen, Phys. Rev. D 22 (1980) 1882. 

[21] V. Acquaviva, N. Bartolo, S. Matarrese and A. Riotto, Nucl. Phys. B 667 (2003) 119; 
H. Noh and J. Hwang, Phys. Rev. D 69 (2004) 104011; K. A. Malik and D. Wands, 
Class. Quantum Grav. 21 (2002) L65; D. H. Lyth, K. A. Malik and M. Sasaki, JCAP 
0505 (2005) 004; F. Vernizzi, Phys. Rev. D 71 (2005) 061301; D. Langlois and F. 
Vernizzi, Phys. Rev. Lett. 95 (2005) 091303; D. H. Lyth and Y. Rodriguez, Phys. 
Rev. D 71 (2005) 123508; K. A. Malik, JCAP 0511 (2005) 005; F. Finelli, G. Marozzi, 
G. P. Vacca, and G. Venturi, Phys. Rev. D 74 (2006) 083522. 



20 



